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We present a many-body description for two-component ultracold bosonic gases when one of the 
species is in the weakly interacting regime and the other is either weakly or strongly interacting. 
In the one-dimensional limit the latter case system is a hybrid in which a Tonks- Girardeau gas is 
immersed in a Bose-Einstein condensate, which is an example of a new class of quantum system 
involving a tunable, superfluid environment. We describe the process of phase separation microscop- 
ically and semiclassically in both situations and show that the quantum correlations are maintained 
in the separated phase. 



Introduction - Impurities immersed in ultracold atomic 
gases have recently emerged as versatile environments for 
studying quantum correlations in highly controllable sys- 
tems [l| . Demonstrations of single fermions or single ions 
embedded in a Bose Einstein condensates (BECs) have 
shown that such systems are experimentally viable Q, 
and are paving the way to study a plethora of new quan- 
tum phenomena, arising from the interactions between 
the impurity and the ultracold environment. Better un- 
derstanding and control of these interactions are already 
leading to new ideas in quantum information theory [3[. 

Precursors to these highly controllable hybrid sys- 
tems have been mesoscopic mixtures of ultracold bosonic 
gases, which consist of either two different atomic species 
or two different hyperfme states of the same species 
They are commonly available in labs worldwide and have 
allowed to investigate mesoscopic quantum dynamics in 
complex systems and to study new exotic states of mat- 
ter. Such systems have a successful microscopic descrip- 
tion based on a two mode model [5] and can be ap- 
proximated semiclassically by using a system of coupled 
Gross-Pitaevskii equations (GPE) [6[. The latter model 
allows to describe the stability of the multi-component 
system [7] and the process of phase separation 0, [![. 

In this work we consider a small number of atoms con- 
fined in an effectively one-dimensional parabolic trap, 
which are either in the weakly or in the strongly cor- 
related regime. Because of the reduced dimensionality, 
the strongly interactingcase corresponds to the Tonks- 
Girardeau (TG) limit [9|, Efl|, which has recently been 
demonstrated experimentally [U|. The same trap con- 
tains a second species of atoms, which we consider to 
be in the weakly correlated regime (the BEC limit), and 
which acts as a tunable environment for the first com- 
ponent. We show that the quantum correlations can be 
tuned through the coupling with between the two species, 
and describe microscopically the phase separation pro- 
cess that drives the immersed species to the edges of the 
BEC. 

A strongly interacting quantum gas in one-dimension 
can be successfully described in different ways. The first 
is to use a model of hard sphere bosons and employ a 



mapping theorem to a non-interacting Fermi gas [9] . This 
permits one to derive an analytical expression for the 
wave function in position space and only relies on the 
knowledge of the solutions of the single particle prob- 
lem H Ell El]- An equal, but numerically expensive, 
description is to use a many-body Hamiltonian and ex- 
pand the field operators into a sufficient number of mo- 
menta [131 ]. 

Here we choose the second approach and expand the 
second quantized field operators for both species in a ba- 
sis of harmonic oscillator functions, including as many 
momenta as needed. This allows us to microscopically 
describe the stationary solutions and the phenomenon of 
phase separation when both species are in the weakly in- 
teracting regime (the BEC-BEC limit) and when one is 
in the strongly interacting regime (the BEC-TG limit). 
While a larger number of momenta is necessary to de- 
scribe the TG gas, only a few momenta are necessary for 
the BEC species, which makes the numerical approach 
possible. We obtain the many body version of the phase 
separation criteria and show that this process occurs as 
a consequence of the excitation of atoms to higher har- 
monics. The phase separated situation can be treated as 
a double well potential for the phase separated species 
and we show that the atoms in the strongly interacting 
species keep their correlations despite being separated by 
the BEC component. 

Since the above approach is limited to small particle 
numbers, we generalise the obtained results to the meso- 
scopic limit by using a semiclassical approach similar to 
the well known coupled GPE systems 0. In the BEC-TG 
limit, however, we show that a GPE coupled to a mean 
field equation with a quintic [14| rather than a cubic non- 
linearity appropriately describes the semiclassical limit, 
while the non-linearity of the cross term depends on the 
value of the coupling between the species. 

Model - We consider a mixture of two ultracold 
bosonic species, confined in a quasi-one dimensional (ID) 
parabolic trap. One component (the 'environment', E) 
is always in the weakly correlated regime, whereas the 
other component (the 'system', S) can be either weakly 
or strongly interacting. The intra-species coupling con- 
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stants are then given by gE,s — —2fr 2 /mE,s a iD, where 
rriE,s are the respective masses of the two components 
and af/ = (-(af S ) 2 /2a E , s )[l - C(a E ,s/af s )]. Here 
cle,s are the s-wave scattering lengths and af_ ,S = 
y / 2h/rriE,s^±^ where u± is the trap frequenc y in the 
radial direction. The constant C is evaluated in [15]. For 
simplicity we assume that the two components are differ- 
ent hyperfine states of the same atomic species and there- 
fore vfiE = vns = M. The coupling constant between 
the two components, gES, is given in a similar fashion 
in terms of the ID inter-species scattering length, af£. 
We expand the field operators in the second-quantized 
description of this system [5| in terms of the eigenfunc- 
tions <j) n {x) of the harmonic potential V(x) = ^Muj 2 x 2 
and use tie,s modes in the expansion for the environment 
and system component, i/je(x, t) = Yln=i a n(t)4>n{x) and 
ifjs(x,t) = Y^n=i bn(t)<f) n (x). The creation and annihila- 
tion operators a\ and a k satisfy the standard (equal time) 
bosonic commutation relations [afc,aj] = Ski, [o^k^i] — 
[a^,a]"] = (and similarly for b\ and b k ), which yields 
the Hamiltonian H = He + Hs + -Hint ? where 

H E =J2 a k a l R kl + \Y1 a \ a \ a m^nVimn (1) 
k,l klmn 

H S =Y J b{kH kl + \Y. b Wi b mb n V k s lmn , (2) 

k,l klmn 

H ™t =- ^2 a \b\brna n V k E i^ n , (3) 



klmn 



with 



H ki = J dx cj)* k (x)H sp (j)i(x), 



(4) 



V klmn = ^T I 
V Mmn =9ES I dx (x)</>* (x)</) m (x)(/) n (x) . (6) 



Here H sp is the single particle Hamiltonian for the har- 
monic oscillator. We then expand the ground state 
= Yli=i Ci ^ i as a sum over a ^ Fock vectors given 

by 



=D 



(^..(^(^...(^-•o, (7) 
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Figure 1. BEC-BEC regime without (upper row) and with 
(lower row) coupling between species. The left column shows 
the average occupations (nt) of the momenta k. Insets are the 
corresponding SPDMs, spanning a region of 6/nm x 6/xra. The 
right column shows the density profiles, with the dashed line 
being the potential, the solid line being the density obtained 
from the SPDM and the solid line with circled markers being 
the one obtained from a semiclassical GPE-GPE simulations. 
The dash-dotted lines in (d) are the approximate Gaussians 
located at the minima of the double well potential Vnl(x) — 
V(x) + gABN B \M^\ 2 l 

BEC-BEC regime - If the interactions gs,s are small 
and gES vanishes it is sufficient to use only a few 
modes to describe both species. Indeed, if the intra- 
species coupling constants vanish as well, the single par- 
ticle density matrix (SPDM), defined as pe(x,x') = 
Y^k,k> (f>k(x)(f)^(x f ) {a\a k >) and similarly for p s (x,x'), is 
simply a Gaussian of width a = ^/h/Muj. For weak 
interactions the atoms still mainly occupy the lowest en- 
ergy eigenfunction ( N E,S w N E,S ) and a Gaussian ap- 
proximation to the eigenstates will be good. Assuming 
spatial overlap between both components, the energy of 
the ground state \I>o can then be found as 



£h°m ={Ne + Ns)^ + -^N E (N E - 1) 



9s 



2V27ra 



N S {N S -1) + 



9es 



2ira 



N E N S , (8) 



with D = (N E l . . . N E JNf\ . . . N^ s l)~i and $ being where we have used that V^o = gE,s/2V^ra. If the 



the vacuum. Here Nf , . . . , N E E (iVf , . . . , N£ s ) are the 
occupation numbers of the tie (ns) modes for the envi- 
ronment (system). The dimension of the Hilbert space is 
n = n E tts with n E ,s = (N E ,s^n E ,s-iy/N E ,sK n E,s- 
1)! where N EjS is the total number of atoms in each 
species. The fast growth of this space for larger particle 
numbers or modes is the biggest challenge to the numer- 
ical approach. 



two species phase separate, one of them, say E, will stay 
in the center of the trap, approximately with a Gaussian 
shape of width ole- A single Gaussian mode then suffices 
in the expansion of the field operator of the environment. 
The system component will approximately assume a pro- 
file formed by two Gaussians of width ^ (see Fig.QJd)), 
which we call (f^ w '. In the harmonic oscillator basis the 
expansion of this state requires a larger number of modes, 
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however, if we use the approximation of <p$ w as two dis- 
placed Gaussians in a potential of frequency uo' , we find 
the energy of this state as 



--Nj 



huo 



9e 



2na E 
9s 



2V2Tra s 



N E (N E - 1) 
Ns(N s - 1). 



(9) 



To determine the point at which phase separation hap- 
pens, we assume that the atoms occupy the same total 
volume before and after phase separation, a E + as = a 
[H, and minimize the energy with respect to a E and 
as- This gives a E = and as = , with 

k = yj ^ ^ ^(Ne~-i) ' Subtracting both energies we ob- 
tain the phase separation criterion 



9esN e N s > ^I^VNe{N e - 1)N S (N S - 1), (10) 

where we have assumed that the change in kinetic energy 
hNs(uo f — uo) is negligible compared to the overall change 
in interaction energies (we justify this assumption below). 
This criterion resembles the semiclassical one for large 
particle numbers, N EiS (N EiS - 1) « N% iS @, Note 
that it predicts phase separation for any nonzero value 
of gES if one of the components consists only of a single 
particle. If we expand uo^ w in terms of the eigenfunc- 
tions of the harmonic potential uo^ w = XI ( ( Po W \ ( t ) k) fa> 
we can see that the process of phase separation coincides 
with the occupation of more and more orbitals in the har- 
monics basis. Consequently, the phase separated species 
will be the one with the higher coupling constant and 
the smaller number of atoms and for simplicity we there- 
fore assume in the following gs > g E and N E ^> Ns. 
This assumption also makes numerical calculations pos- 
sible, since the shape of the environment component does 
not change substantially after phase separation (a E w a) 
and it can correspondingly be described only considering 
a single mode. 

Let us consider an example using atoms of the mass 
of 8T Rb in a trap of frequencies uo = 2tt • 400Hz and 
uo± = lOOuo. We (exemplary) choose the scattering 
lengths to be a E = ao/10 and as = 10ao, where ao is 
the Bohr radius. The environment consists of N E = 300 
particles and the immersed component of Ns = 4, giving 
g E = 10~ 9 ftw and gs = 1.01 x 10~ 7 ftw. In all calcula- 
tions we use ns = 14 and n E = 1 (as justified above). 
In the upper row in Fig. [T] we show the situation where 
the different species do not interact (g E s = 0), while 
in the bottom row g E s = which is deep inside the 
phase separated regime. As expected, for g E s = the 
SPDM of both species is Gaussian, as shown in the inset 
of Fig.[TJa) for the system component. The energy of the 
sample in this case is given by eq. (|EJ) as £; hom = 187.0 huo. 
Correspondingly, the average occupation for the system 
component (rik) = (b\bk), represented in Fig.QJa), shows 



that mainly one momentum component is occupied, be- 
cause the coupling constant is very small. For bigger gs 
higher lying momenta will start gaining occupation, how- 
ever the system will still form a BEC as the occupation of 
the lowest momentum component is of the order of Ns. 
As long as the sample is in this limit no phase separation 
occurs. Fig. []Jb) shows the density of the system com- 
ponent obtained from the SPDM (solid line) coinciding 
with the semiclassical calculation using two coupled GPE 
(solid line with circled markers). 

If 9es — (lower row in Fig. [1]) the energy for 
a phase separated state is much smaller than that for 
a mixed phase, E inh = 191.0 huo < E hom = 230.8 huo 
and also corresponds approximately to the one calcu- 
lated numerically of E = 195.0 huo. Accordingly, the 
SPDM in this situation (see inset of Fig. [He)) shows 
that phase separation has occurred and an occupation of 
higher momentum modes is found (see Fig. Die)). The 
densities calculated from the SPDM and obtained from 
the numerical solution of the coupled GPE shows again 
good agreement (see Fig. 0Jd)). The Gaussian functions 
used to calculate the energy given by eq. (j9j) are located 
at the minima of the double well potential Vnl(x) = 
V(x) + gEsNs\fa(x)\ 2 \, which are given by d = ±av7?, 
with f3 = In {2g ES Ns/a 3 M^/7Tuo 2 ). The local trapping 
frequency can then be approximated as uo' = 
leading to a width of the Gaussians of a' = ^/h/Muo' . 
We can now estimate the increase in the kinetic energy 
due to phase separation as hNsuo(l — a/2/3), which has 
to be compared to the change in the interaction ener- 
gies (g ES N E N s and y/g^y/N E (NE - l)N s (N s - 1)). 
Since the latter is quadratic in the number of particles 
it is generally much larger except for systems with very 
small coupling constants and number of atoms. 

BEC- TG regime - We now consider the case where the 
system component is in the To nks- Girardeau regime. In 
Fig. [2]we show the same quantities as before, but now for 
N E = 40, a E = ao and as = 500ao, which gives a Lieb- 



Liniger parameter of 7 = 2gsML/h 2 Ns = 4.5 [16|, where 
L is the size of the cloud. This corresponds to the system 
component being in the To nks- Girardeau regime and we 
find the interaction coefficients to be g E = 10~ 8 hj and 
gs = 9.6 x 10 _5 fkj. As expected, for no inter-species in- 
teraction the SPDM and the momentum distribution for 
the system component resemble that of a TG gas (Fig. [2] 
(a)). The density profiles calculated from the microscopic 
and the semiclassical approach (see below and Fig.[2](b)) 
show again good agreement. Though the phase sepa- 
ration criterion of eq. ([TP]) is only approximately valid 
in this regime, due to the deviations from the Gaussian 
shapes for the individual components, it is still useful 
away from the precise border transition point. The lower 
row of Fig. [2] shows the situation for g E s — 0.05^5, where 
phase separation is clearly visible in the SPDM and the 
momentum distribution shows a shift of the momenta 
due to the density resembling a higher excited state. 
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Figure 2. BEC-TG regime without (upper row) and with 
(lower row) coupling between species. Panel lay-out as in 
Fig. [T] The semiclassical profiles are obtained using the mean 
field equations described in the text. 

To compare the above results with a mean-field 
model in the semiclassical limit, we chose to model the 
strongly correlated system component using a quintic 
non-linearity in the field equation 



14j . While one has 



to be aware that this approach for a TG gas does not de- 
scribe the coherence properties correctly [17j , it is known 
to give a good approximation to the exact density in the 
single component case. This equation is then coupled to 
a GPE for the BEC component, giving 



itys = ~ i^'s + [V(x) +g s \ips\ 4 + 9es\^e\ p ] ^, 



s 

(lla) 



h 2 ... 



ihipE = - -^Pe + [ v ( x ) + 9e\iPe\ 2 + 9es\iPs\ p ] 4>E, 

(lib) 



where gs 



2M 



. If the interactions between both com- 



ponents are small, i.e. j' = 2gEsML/h 2 Ns < 1, the ex- 
ponent on the non-linear coupling term is given by p = 2 
and in the opposite limit by p = 4. For the parameters 
used in Fig. [2] we find 7' = 0.22 and the numerical so- 
lution of eqs. (|TTJ) show good agreement with the micro- 
scopically calculated density (see Fig.EJc)). Calculations 
for increased values of gss, which require p = 4 also show 
good agreement, but are not shown here. 

Let us finally discuss the pair correlation functions 
g( 2 \x,x f ) = (&(x)&(x')4/(x f )^(x)}, with 4f(x) being 
the field operator, using again the microscopic model. 
These correlations give the probability of finding an atom 
at position x once another atom has been measured at x' . 
The BEC-BEC case with no interaction between species 
(corresponding to the upper row of Fig. [T]) is shown in 




x' x' 

Figure 3. Pair correlation functions g^ 2 \x,x / ) spanning an 
area of 6/ira x 6/ira. The top row shows the BEC-BEC cases 
using the same parameters as in Fig. [T] for the situation (a) 
without and (b) with coupling between species. The lower 
row shows the BEC-GP regime analogous to the situation in 
Fig. [2] for the situation (c) without and (d) with coupling 
between the species. 

Fig. [3fa) and the expected Gaussian profile is obtained. 
In the phase separation regime (see Fig. Efb)) we find 
that the two approximately Gaussian parts of the density 
profile of the system (shown in the lower row in Fig. [I]) are 
highly correlated . For the BEC-TG case with no inter- 
species interaction (corresponding to the upper row in 
Fig. [2]), as shown in Fig. [3jc), we recover the well-known 
hard-sphere behaviour that no two atoms can be found 
at the same point in space, which persists in the phase 
separated limit for the individual peaks (see Fig. Efd)). 
The correlations between the two peaks are also still visi- 
ble, however they show a more complicated structure due 
to the inter-species interaction. 

In conclusion, we have presented a microscopic and a 
semiclassical model to describe a two-component Bose 
gas at ultralow temperatures in the one-dimensional 
limit. Since we allow for one component to be in any 
correlation regime, our work generalises the known mod- 
els for interpenetrating Bose-Einstein condensates to in- 
clude a significantly larger group of systems, which are 
currently about the become experimentally available. 

From our microscopic model we have, as a first ex- 
ample, derived a criterion for phase separation, which 
extends the well-known mean field result to the meso- 
scopic limit. We have also presented a semiclassical de- 
scription of a mixture of ultracold bosons when one of 
them is highly interacting, and found that the nonlinear 
coupling term between the species depends on the inter- 
action strength. Finally, we have shown that even when 
the atoms are phase separated, they are still strongly cor- 
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related among themselves. Interesting questions that can 
be approached with this model include the study of quan- 
tum correlations between a system and an environment 
in fundamental settings as well as using one component 
(one matter wave) to engineer the state of the second 
matter wave. 

This project was support by Science Foundation Ire- 
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